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We propose a Kaluza- Klein cosmology by reversing the usual scenario: instead of starting 
with a flat 4-1- N dimensional universe in which N of the dimensions curl up into a compact 
manifold, we start with a compact 3+ N dimensional manifold in which 3 of the dimensions are 
allowed to peel off and expand into the known universe. 

We reverse the usual "spontaneous compactification" scenario and begin with a closed manifold 
M 3+ " which undergoes "spontaneous fissioning" into a product manifold M 3 xM N . Remarkably, 
the 3-dimensional universe M 3 can undergo a rapid de Sitter expansion large enough to solve the 
horizon and flatness problem. We call this "topological inflation", which we propose as an alternative 
to the usual GUT inflation. The inflationary phase automatically terminates into a big bang phase. 



1. Introduction 

The Kaluza- Klein theory [ 1 ] is one of the most promising and attractive candidates 
for a unified field theory which includes both the gravitational and particle interac- 
tions. By embedding Einstein's theory into a higher 4+ N dimensional manifold, 
the Yang- Mills interactions emerge naturally when we decompose the components 
of the higher dimensional metric tensor. Although the Klein- Kaluza picture is 
appealing from the gauge theory point of view, the application of Kaluza-Klein 
theory to cosmology remains totally obscure. 

In the usual Kaluza-Klein scenario, we assume that the universe was originally 
defined on a flat 4+N dimensional manifold, in which N of the dimensions 
underwent spontaneous compactification [2], i.e. N dimensions suddenly curled up 
into a compact manifold the size of the Planck length, leaving a flat four-dimensional 
riemannian universe. 

In this paper, we apply Kaluza- Klain theory to cosmology by reversing this 
scenario. Instead of N dimensions curling up into a compact manifold, we assume 
that the original universe was defined on a 3 + N dimensional compact manifold 
M 3+N on the scale of the Planck length which suddenly "fissioned" into a product 
manifold M 3 x M N . The M 3 universe then exploded outward into the known universe, 
while the M N universe collapsed. 
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The usual Klein- Kaluza picture assumes spontaneous compactification, where a 
non-compact manifold suddenly becomes balled up into a compact manifold. We, 
however, assume the opposite— 'spontaneous fissioning [3], i.e. the original compact 
manifold "fissioned" into two smaller compact manifolds, one of which then 
exploded rapidly into the known universe while the other compact manifold sudden 1 v 
collapsed down to the Planck length. 

Remarkably, we find that the rapid expansion of the 3-dimensional universe is 
in an exponential de Sitter phase. With a suitable choice of parameters, we can 
achieve sufficient exponential expansion ("inflation" [4]) in order to solve the 
long-standing horizon and flatness problems. 

We propose this "topological inflation" as a solution to the horizon and flatness 
problems. Because the usual GUT inflation picture has certain problems with fitting 
or fine tuning the parameters [5], this "topological inflation" formalism is an 
attractive alternative to the usual formalism. 

We find that the fissioning takes place in three stages: 

(I) In the first stage, the original compact manifold M 3+N undergoes a de Sitter 
expansion in 3 + N dimensional space. This manifold then suddenly fissions into u 
product manifold of M 3 x M N initiated, for example, by the Freund-Rubin ansatz 
[6]. 

(II) In the second stage, the three dimensional manifold undergoes a rapid 
exponential expansion, while the N-dimensional manifold collapses down to the 
Planck length. This period of topological inflation can be made sufficiently large to 
solve the horizon and flatness problems. 

(III) In the last phase, the de Sitter expansion automatically shuts itself off and 
smoothly continues over into a big bang phase (sec fig. 1). 

The attractive features of our model are as follows: 

(i) We have checked explicitly on computer that the M 3 universe undergoes an 
inflation phase which eventually cuts itself off and goes over into a big bang phase. 
A suitable choice of parameters can create a sufficient amount of inflation to generate 
the desired e 65 factor. 

(ii) Although we have to fix the value of the primordial cosmological constant 
before spontaneous fissioning takes place, we find that our formalism predicts an 
extremely small but non-vanishing value of the effective cosmological constant when 
viewed from the present universe. This may help to shed some light on the cosmologi- 
cal constant problem. 

(iii) Supersymmetry can be probably included in this model if we use the 0(8) 
model to generate spontaneous fissioning. If we calculate quantum corrections to 
an 0(8) ansatz which preserves supersymmetry, however, first-loop corrections 
vanish between bosonic and fermionic contributions [7]. Therefore, we must calcu- 
late the quantum corrections around a supersymmetry-breaking ansatz. Also, because 
our formalism generates a series of terms which contribute to the cosmological 
constant, this may shed some light on the problem of how to convert the supersym- 
metric anti-de Sitter phase into a de Sitter one [8]. 
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Fi« 1 (a) In re R ion I the universe is a closed manifold in 3+ N dimensions. At point A the universe 
h tfiss otd inTo'Tc.; Ld 3-dimensiona. manifold and a N-dimensiona. manifold. The 
Z£SZ£&m a rapid de Sitter expansion, while the 
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(iv) Although we performed our calculation on closed spheres S N , we could also 
have used other homogeneous manifolds such as CP". Using complex manifolds, 
we may be able to incorporate GUTs into our formalism, although it is independent 
of GUTs. 

(v) Although we use a zeta functional regularization to define our theory, we find 
that our equilibrium results agree remarkably well with the results of Candelas and 
Weinberg [9]. 

We should make perfectly clear in our paper that the precise mechanism by which 
we break the initial symmetry of the universe, whether by Freund-Rubin type 
mechanisms or otherwise, is left completely unsolved. This is the fundamental 
weakness of this paper. As a result, it is difficult to argue precisely how much 
inflation actually takes place in this model, because in principle any initial values 
of the radii are possible after fissioning. We have, in a sense, simply pushed back 
our understanding of the fine tuning problem from the parameters found in the 
GUT model to the parameters of the original dynamical breaking of the Kaluza- 
Klein universe. 

As a result, we should be careful to state that we certainly have not solved the 
fine tuning problem. What we are claiming, however, is that we are postulating a 
new mechanism which (in principle) may solve the fine tuning problem but which 
(in practice) depends crucially on how the original symmetry of the Kaluza- Klein 
universe was broken, which nobody knows. We simply wish to correct the conven- 
tional wisdom, however, that inflation takes place solely via GUT mechanism or 
even super-GUT mechanisms. However, we also wish to state that it must be 
physically possible to experimentally justify the claim that our Universe once went 
through a Kaluza-Klein universe, and that one of these physically relevant 
phenomena might be Kaluza- Klain-induced inflation. 



2. The fissioning universe 

Our starting point is the usual Einstein theory defined on a closed 3 + N 
dimensional sphere with a cosmological constant: 

RAB-jgAi,(R + A) = -STrGT AB , A, B = l,2 4+N. (2.1) 

Let us assume that we can spontaneously break the original manifold of the theory 
down to S 3 xS N . We assume that we can now split the metric tensor of the theory 
into the following: 

gAB = (-l,gy,g a/3 ), g,„=0, 

U=l,2,3, a, /3 = 1,2 N. (2.2) 

Now let us assume a Robertson- Walker ansatz defined over this product manifold. 
Let us assume that the 3 (N)-dimensional sphere has a radius given by p(p). We 
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easily find that: 

Roo=3p'/p + ^P/P» 

« y = *, - g y (p/p + 2(P/P) 2 + N(p/ P )(p/P)) , 

R oP = ^-^(p7p+(N-i)(p/p) 2 +3(p7p)(p/p)). 

r = £ - 6p/ p - 2 NpV p - 6(p7 p) 2 

-N(N-l)(p/p) 2 -6N(p/p)(p/p)- (23) 

Inserting these identities into Einstein's equations, we find: 

3(p/p) 2 + 5N(N-l)(p/p-) 2 + 3N(p/p)(p/p) + 3/p 2 

2p7p+(p7p) 2 +Np7p+5N(n-i)(p/p-) 2 +Vp 2 _ 

+iN(N-D/p 2 + 2N(p/p)(p/p)-M = -5^^g v , 
3p7p + 3(p/p) 2 + (N-1)p7p + KN-1)(N-2)(p/p) 2 + 3(N-1)(p/p)(p/p) 

+ i(N-l)(N-2)/p 2 -M 

_ 2 8tt = (2.4) 
+ 3/p "n opg ' 

Now let us calculate the contribution to the energy-momentum tensor. By 

scalar field, which is now defined over the product manifold. 

|^ c ^> = der ,/2 [-3 1 +V 2 +m 2 (t)]. ( 2 - 5 ) 
(The factor m(t? is a function of the bare mass as well as certain time derivatives 

^should note that (2.5) and our basic conclusions are based on 

*h»t the "fissionine" process has already termmated and the radius 01 me 
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Planck times have already lapsed. This is because of two reasons. First, quantum 
effects will dominate the early era when the two radii are of equal magnitude. Thus, 
our entire mathematical framework breaks down in this early period near t = 0. 
Second, we wish to avoid the problem of time-dependent adiabatic effects. At later 
inflationary times, when the three-dimensional radius is much larger than the 
TV-dimensional radius, time-dependent effects are explicitly included in m(t) in 
(2.5). We have calculated such effects, and find that they can be neglected if the 
radius of the three-dimensional universe is much larger than that of the TV- 
dimensional universe. But this is precisely the physical domain which we are 
studymg. In other words, the time-dependent terms in m(t) in (2.5) have been 
calculated and can be shown to be small in the region we are studying, when the 
radius of the three-sphere and the logarithmic time derivative of the TV-sphere are 
both larger than the radius of the TV-sphere. Thus, we can drop such terms in this 
paper. (See the last paper in ref. [9].) 

We now exponentiate this determinant in order to calculate its contribution to 
the energy-momentum tensor. 

The determinant over the product manifold can be computed as the product over 
the various eigenvalues of the laplacian. On an TV sphere, the ith eigenvalue of the 
laplacian is given by /(i + TV-1) with a degeneracy of (2i+/V-l)x 
(«'+TV-2)!/(/V-l)!/!. 

In order to put in the temperature dependence into our formalism, we assume 
that the scalar field is defined over a periodic time interval with a period given by p : 

fUO-Ze^Mx). (2.6) 

We are making the assumption that the scalar field is a periodic function defined 
over a finite time interval because it is in thermal equilibrium in a background 
Robertson- Walker metric. We assume that the background metric, however, is not 
in equilibrium, i.e. it expands in a normal background de Sitter phase without 
periodicity. (Introducing temperature into our formalism via periodic functions 
defined in a non-periodic background gravitational field is also used in the usual 
GUT formalism.) 

Crucial to this assumption is the condition that quantum gravitational corrections, 
which are of the order of the Planck length, can be neglected. We can implement 
this assumption if there are, say, several hundred or thousand matter fields coupled 
to gravity. The dynamics of fissioning and inflation then take place in an energy 
region several orders of magnitude larger than the Planck length, so we can con- 
sistently neglect quantum gravitational effects. Our approximation, however, breaks 
down in the high-temperature domain, where we can no longer neglect quantum 
gravitational effects. Therefore, we have focused our discussion on low-temperature 
effects rather than high-temperature effects, where our approximation may actually 
break down. 
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Putting everything together, we find that we can express the energy-momentum 
in terms of a function V given by: 

V = — 2 T (i + \) 2 (2j + N-\)(j+N-2)\/(N-\)\j\ 

Jk -co 

xln M^ + W±^ + i£l. (2.7) 

l P P P J 

(The prime denotes the deletion of the i = j = fc = 0 term in the summation.) 
We find that the components of the energy-momentum tensor can be written as 
derivatives on the function V: 

Toog m, --(p 3 p N i3)-^^V, 

T,g^"(p 3 p N ^rVd P v, 

T^-Vp^aiV, (2.8) 

3. Zeta function regularization 

The function V defined earlier formally diverges when we take the sum. In order 
to regularize the function V, we will find that the zeta functional method is the most 
convenient of several regularization schemes. 

Let us define the following function: 

£' (i+l) 2 {2j+N-l)(j+N-2)\/(N-l)lj\ 

k = -oo 

IP p p J 

In terms of this function, we can re-express our function V, and hence all 
components of the energy-momentum tensor: 

v = _J l(r (0) + f(0)ln/i 2 ), (3.2) 

4 77 

where /i is an arbitrary mass scale. 

At this point, the function still diverges in the region of interest. Our next step 
is to re-express the function in terms of the generalized Epstein zeta function [10], 
which has a known analytic continuation. Let (g,, g 2 • * * g»), (*i, h 2 • • * M and 
(m,, m 2 • • • m n ) represent sequences of n numbers. Let c i} be a symmetric n x n 
matrix with an inverse. 
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Let us define the nth rank Epstein zeta function as: 

Z (h h ... h) is)s * ••• I e 2 -^.'".V l/2 ((m + g)), 

*((*))"■£ CyX,x 7 . (3.3) 

It is now straightforward to re-express the function f totally in terms of derivatives 
of Epstein zeta functions of third, second, and first rank: (for convenience, we set 
N = 3, we rescale p -* 2irp and we let the prime indicate that the lowest term in the 
summation for m, n 9 k is deleted; if k goes from minus infinity to positive infinity, 
then the prime indicates that k = 0 is deleted) 

m,n = l \p p & J 

fc = -oo 

k = -ao 

4 —4 ~* 



where 



, _m 2 » 2 it 2 »» fc 2 

m 2 k 2 ,,11 
fmk^-^+-i=9„p, b 2 = -^+j 2 -m 2 {t), (3.5) 

A 0£±£L_ 6 2r f r ,- r .l T(5 + r) , f /0 0 o\ 

4 Wi+r)* i, r "*-4r(,)ra+r)' rlo o o/J 2s+2r) 

~ Z (o o)^ (2s+2r)_Z (o o), ( 2s + 2 '') + 2 / 32l+2, ^ + 2r)-4(fc 2 )— }, 



(3.6) 



V V 9 dp 2 ) r(*-i)r(r+i)-r\o o ' 

~[ Z (o o)„ (2* + 2r-2) + (p^p)]+2/3 2 ' +2r - 2 i:(25 + 2r-2) 

-4(fr 2 )— ( 3. 7) 
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' 2 ((m + g)h 

(3.3) 

ms of derivatives 
venience, we set 
west term in the 
positive infinity, 



■b 2 y* +i 

C„ (3.4) 



(3.5) 



(2s+2r) 



(3.6) 

! + 2r-2) 



(3.7) 



i P y a a r(s +r -2) ( Jo o o\ 

Cr 4(s-2)(s-l)ap 2 <9p 2 r(s-2)r(l + r) [ \0 0 0/„ v 

~[ z (o o) ( 2i+2 '- 4 ) + (p^p)] +2 ^ I+2r " 4 ^ 2s+2r - 4) - 4(b2) " I_r+ i- 

(3.8) 

Notice that when s = 0, the function f diverges. In order to regularize this function, 
we now make use of the following relations, which allows us to analytically continue 
an nth rank zeta function defined at the point s to a point n-s: 



= 4 1/2 exp [-27TI 1 g t h^ir~ in - x)/2 rti(n-s)) 

I h y h 2 • • • h n \ 
xZI (n-s), 



(3.9) 



^-det(c 0 ), 



3% 



(3.10) 



At this point, we will still need one more formula, which is the value of the zeta 
function at s = 3 and s = 2. In the appendix, we will prove that the zeta function of 
fixed rank at a given point can always be re-written in terms of the sums of zeta 
functions of lower rank. Because zeta functions of first rank are the usual riemannian 
zeta functions, this means that an Epstein zeta function of arbitrary rank can 
eventually be reduced to sums over riemannian zeta functions. In the appendix, we 
will show that: 



»-2 \o o/, 



(s) 



2ir it 2iry a 
4 1/2 (s-2) 3a A 1 ' 2 4A 



~ I ln(l-a 2 ")(l-a!"), 



(3.11) 



lim ZK ^ *J 
*-3 \0 0 0/. 



(s) 



1* , 2f(3) | 2jr_ | 4ir(y-l) 47rln)3 8tt - , n (1 _ e -2^»/ P) 

s-3) p° 3p 2 p ppfi ppj8 pp/3 B =! 



4ir 

pp/3(s-3) p" 
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\fl 00 J 1 1 1 

— Z | ii" 3 |^ln(l-e- 4 ^)+^In(l-e- 4 ^)| 
4y/2 00 f 1 1 1 

P «-i lp p 2 J 

«v/2 00 

P m,*=l 

2v/2 00 

— ^I i (p a « a +^r va InO-c- J, -) i (3.12) 
¥> 2 ■ am? + 2bm, m 2 +cm\ , 4 = ac - b 2 , 

B mk = — (p 2 m 2 +p 2 k 2 ), ym 0.57721566.... (3.13) 

We now have all the tools necessary to re-express the function V entirely in terms 
of simple summations and known analytic functions. The final calculation is straight- 
forward but quite tedious. We will present the final result in the appendix. 

We will find, however, that in the regions of interest, almost all the terms drop 
out and the potential function reduces to a remarkably simple result. 



3. Topological inflation 

At this point, we have reduced the function V to a series of known functions and 
easily computable summations. 

Although the problem of computing V is now formally solved, we will find it 
convenient to make certain approximations which will allow us to compute the 
cosmological behavior of our equations. 

We assume self-consistency in region II, that the radius of the 3-dimensional 
universe is much larger than the N-dimensional universe. We also assume that the 
mass terms are small so we can drop the mass term m(t) in our potential (we will 
return to this point later). In this approximation, we can reduce the expression for 
the zeta function to the following form: 



p-o \o o o/. 1 ; 



=2pV(^)+2>^r((5-i)/5)r-H^)p/3-'j:(5-i)+ pm ^" 2 y j) p-'-V(5-2) 

+5 n£r (Ww)<5 " 2)/2 ' K< ^ 2,/2(2,r/3/cM/p " ) • w 
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limit of large p: . 

A, .hi, point, w. can « begin the search for oosmoiogic, solurions ,0 the 

equations of motion. magnitude smaller than the 

/» region I, where the time scale is ma y ( , deSitter expansion. 

Planck length, we find that the mamfold M h to .civ,,, 

v = c 3 p 3 /3/p 4 -^ 2 (p//3) 3 . 

C 3 = 7.5688---xl0- 5 . (4 ' 3) 
This compares f.vorah.v wUh ,he C^^:^ 

Given this simple expression for V, we can now v 
of motion into the final form: ^ 

3x 2 + 3/p 2 -^ + [5/+7+9xy + 3/p 2 ] = 24.Gp^ 2 ^ 

3x + 6x 2 + [3/ + 2y + 3/p 2 ]-U + l/P 2 = ^Gpi(-4C3/p 7 ), 

XS p/p, y^plP- 
,„ region .,, * us assume that and T*- -« 'actor, in the 

brackets vanish in (4.4). ^ analoev with the usual GUT 

A solution for these equations ex.sts ^^Tom ^Z equations generated 
picture, in many ways, these equations mi m.c the GUT mfla on eq ^ 

by Higgs particles. In some sense, *e , ^hey a ^simple newtonian 

resemble the equations of motion for the H.gg part^^ 



like terms.) 
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Using the newtonian-GUT picture as a guide, we can find the following solution 
for the equations of motion in region II, where the M 3 radius undergoes a de Sitter 
expansion while the M N radius is fixed: 

P = (1.46033- ■ -)p 0 , 

P ~exp [0.515305- • (r/p 0 )], 

y = y = 0, )8-' = 0, 

G = Gp 3 0 , p 0 . (0.047109 • • • )G m . (4.5) 

Notice that this de Sitter expansion eventually terminates on its own. When the 
radius of the ^-dimensional universe begins to move away from its value of 
(1.40633 • • • )p 0 , it gradually approaches its equilibrium value at p 0 . When the 
radius of the JV-dimensional universe finally reaches its equilibrium value, then the 
de Sitter expansion terminates. 

In region III, the three-dimensional universe undergoes a simple big bang driven 
by the matter content of the energy-momentum tensor. The JV-dimensional universe 
has now settled down to its equilibrium value of p 0 while the three-dimensional 
universe undergoes a simple power law expansion (big bang) 



P~t m , 



P = Po, 



^ - TPo . 

p~t 1/2 . 



(4.6) 



We have checked explicitly on computer that the inflation phase can be made to 
extend for arbitrarily long periods of time and that it eventually terminates into a 
big bang phase. 

Notice that the value of the primordial cosmological constant used in region I is 
now determined by the equilibrium equations of region III. In order for the inflation 
to terminate smoothly into the big bang phase of region III, it is necessary to fine 
tune the cosmological constant at the beginning of time. Our picture, therefore 
does not give a solution to the cosmological constant problem. Once the cosmological 
constant is fixed at the beginning of time, however, our scenario explains why the 
effective cosmological constant at the present time is so small. Because the Planck 
length is the only dimensionful parameter in the theory and because considerably 
many Planck lengths have passed since the original fissioning, the effective cosmo- 
logical constant found today must be astronomically small. Our picture then 
actually makes a prediction as to the present value of the cosmological constant' 
Unfortunately, the precise value of the effective cosmological constant at the present 
time is too small to be reliably calculated. 

^The next question we ask is how many parameters are necessary to achieve the 
e factor necessary to explain the horizon and flatness problems. In the usual GUT 
picture, there are certain difficulties in fine tuning the parameters to achieve the 
proper e expansion factor. 
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POTENTIAL 




Fig 2 As in the case with GUT inflation, the evolution of the "topological inflation" phase depends 
on the time it takes for a ball to roll down a potential curve. Because we have the freedom to choose 
the values of the initial parameters, such as the mass of the scalar particles, we can flatten th.s curve 
and vary the time period over which inflation takes place. (Without a proper choice of initial parameters 
however, the ball would have to be placed directly on top of the hill, at pointB, «c .insure ha, ^we , e 
enough inflation. This can conceivably reintroduce the fine tuning problem.) The shape of the curve is 
quite dependent on a variety of factors, such as the mass. Unfortunately, we do not 
mechanism which determines how the universe fissions, and therefore do not have the initial parameters 
which will determine how much inflation takes place. 

If we examine the equations for the radius of the N-dimensional universe, we 
notice that they resemble the equations for a newtonian particle in a potential well, 
similar to the situation in GUT inflation scenarios with Higgs particles. The inflation 
takes place when the radius of the N-dimensional universe is placed close to its 
unstable equilibrium value of (1.40633 • • • )/> 0 , which lies on a hill-top in potential 
space. If we take the mass to zero and set all friction-like terms to zero, then we 
have a simple potential curve as in fig. 2. 

Inflation takes place while the particle is on the top of the hill as well as while 
it rolls down the hill. Inflation terminates when the particle finally asymptotically 
oscillates around the very bottom of the hill. If there are no mass terms m(f), then 
we must precisely fine tune the position of the particle directly on top of the hill, 
which is undesirable. 
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In reality, of course, there are several parameters to play with. First of all, the 
mass term for the scalar field changes the shape of the original potential curve by 
flattening out the hill top and therefore increases the time interval in which inflation 
takes place. 

For small mass, the net effect is to lengthen the time of inflation. The shape of 
the potential is now a function of a continuous parameter m. For large masses, 
however, it may actually be possible to significantly alter the shape of the potential 
curve, including the addition of more than one local minimum and maximum. We 
will explain this in more detail in another paper. 

Second, we have the time dependent terms coming from the original evaluation 
of the 34- JV dimensional laplacian. These friction terms, although they vanish in 
the region of interest (at the top of a local maximum) contribute when the particle 
rolls down the hill. 

Third, we can still adjust the parameters coming in from the original fissioning 
process, i.e. the initial value of the radii and its derivatives when fissioning took place. 

5. Conclusion 

In this paper, we are proposing an alternative to both the usual Kaluza- Klein 
and the standard GUT inflation pictures. When we reverse the Kaluza- Klein 
scenario, and allow three dimensions to peel off from a 3 + AT dimensional compact 
manifold, we find, remarkably enough, that the three-dimensional manifold under- 
goes a de Sitter-like expansion, while the other N-dimensional manifold collapses 
down to the Planck length. By varying certain parameters, we have enough freedom 
to achieve an arbitrary amount of inflation, enough to explain the horizon and 
flatness problems. In a later paper, we will elaborate on precisely how these 
parameters affect the amount of inflation. 

The application of the Kaluza-Klein formalism to cosmology, therefore, leads to 
a consistent picture compatible with the known features of the early universe. 

Our method stands in contrast to other time-dependent formalisms involving 
Kaluza-Klein theory. Our method uses quantum one-loop corrections to generate 
the effective cosmological term and to drive the inflation. Alternatively, one might 
conceivably use external factors, such as higher derivative terms, the equation of 
state for a classical gas, or even an external monopole [11-13] to generalize the 
static Kaluza-Klein picture. We find, however, that these terms unnecessarily compli- 
cate the model, and not all of them are compatible with inflation. Single-loop 
corrections, in any event, are an essential feature of any Kaluza-Klein model 
regardless of external features one may add, and single-loop corrections alone are 
sufficient to give satisfactory cosmological behavior. Single-loop corrections generate 
acceptable cosmological behavior by themselves without adding external features. 

Again, we stress that we certainly have not solved the fine tuning problem. In 
some sense, what we are doing is replacing our ignorance of where inflation comes 
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from with our ignorance of how the original universe was spontaneously or dynami- 
cally broken. Before we can realistically calculate exactly how much inflation our 
scenario produces, we must know the initial radii of the two universes after fissioning. 
This, in turn, implies that we must wait until we understand symmetry breaking 
schemes before we can say with confidence precisely how much inflation this model 
produces. In this sense, we are replacing one set of unknown parameters with another. 



Note added 

After this paper was written, we found that M. Yoshimura [14] also investigated 
the time-dependent behavior of Kaluza- Klein theory when single-loop corrections 
are included. He suggests that the mechanism of oscillating Newton's constant may 
help to understand the horizon and flatness problems in the high-temperature limit. 
Our paper analyzes the behavior at relatively low temperatures, where we find 
inflation. (We avoided an elaborate discussion of the high-temperature limit in our 
work because the fundamental approximation in our work, that quantum effects 
from gravity itself can be dropped, begins to break down in the high-temperature 
limit) 



Appendix A 

In this appendix, we analyze the decomposition of the V potential in terms of 
zeta functions. Using the identities presented earlier in the paper, we can reduce 
the potential function V to the following form: 

V- I A r + B r +C r , (A.1) 

where 

A n _l£(l + l) +ln{ j2rt b >) + I Ind-e-^Kl-e- 2 ^), (A.2) 
12 \p pi »-i 

/2 = ^+T?-*M2^)-I (ln(l-e-^"") 

47rp 12p n=l 

I 62<e*— I pp^ 2 ln(l-e- B -), (A3) 

7T m ,k=l 47r »hk~l 
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A, = ^(° ° J) ( 2 )-(^ Z (o o). (2) + (P-P)) + A^V + lni, 2 ,(A.4) 

1 CO / +00 \ 

A r (r>\)=-b ir r fz k +U r /m£+(p«*p) 

k 00 

"W JL/ m/m * + ^ 3 "2. 3 J./- 4*vVr (4) 

+ (p**p)}+21nfr 2 , (A.6) 

-2 £ Ind-e-W'n-irPpb 2 I /i(e Wp -l)- 1 + (p«p)] 

- £(0) In 0 2 - 2H0) - 2 In b 2 - 2 , (A.7) 

+Vp 3 /3 I /m^l+lnfe 2 , (A.10) 

m t n,k = -cxD ) 

-\{Wflb 2 - P m I' p 4 m 2 f-l' k 2 -&pp 3 f}b 2 -ppl3) 1' p*n 2 rt£ 



m,n,fc*= -00 J 
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Jell Y f~2 1P1£ y m 2 f ~3 

+^a4)p 3 P(p 4 +r 4 )-^f+(p^p))-2lnb 2 -l, (AM) 

(-4 4 £ + ^)H t 1^ , <' I, "-'»<'-«- 4 -"> 

-^pp 2 /3 i »-V^" /<s -i) _, -^pp^ I n-V ,r '" ,/^ -l)" , 

077- n = i 47T * = ] 

V2pp 4 /3 - _ 3 ri67rV _e^^_ 4ira/p 2 ~| 
32w 2 ~," L P 3 (e 4 "" n//i -l) 2 e^'-lj 

P 2 W S , e 4 "^ ,^ ? ..-J*™ e 4 ^ 

V2 .t, (e 4 ^-l) 2 8tt P P * " L p 2 (e 4 ^-l) 2 

(16ff 2 pn 2 /p 2 )e 4,rp "^ 32ff 2 pw e girpn/ n - a » m 2 e fl -> 

v^pVft • , [ 16ff 2 pw 2 e B "-> 32pm 2 e 2B ~* 1 
4tt B ,.i/"" E Lp 4 (e B --l) 2 p 3 (e B ~'-l) 3 J 

I07T m,k = l 4 77 m,/c-l 



16tt 

4 . 4 

16*- 2 



„"=4 fl S f 12ttpw 2 , ,_3 /2 167r 2 pm 2 e B -* 4wpm 2 ] ) 
P 9 P Jfe., lp 2 /^(e B — p 3 (e B — I) 2 "p" 2 / 2 ,(e B — 1)1/ 



+§ln(2*j8)-^U E lnd-e" 2 ^) 
12p „-i 

-fc 2 [v^/p 2 | ] «(e 2 ' ran//s -l)- 1 +^7ri8p + (p^p)] 

-|ln^ 2 -^ln2ir + ln/3 2 +|, (A.12) 
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Appendix B 

In this appendix, we will show how to express an Epstein zeta function of third 
rank in terms of zeta functions of lower rank. This method can be generalized to 
zeta functions of arbitrary rank. 

The starting point is the following identity, which we shall use throughout our 
calculation: 

I e - w = 4= I e-™ 2/2 . (B.l) 

m = -oo \ Z m ^-ao 

We shall also use the following expression for the Epstein zeta function: 

* :::■ j) w 

= dzz l/2 -' 1 • • • I exp -irz<p((g+m))+2m £ m,h, 1 . (B.2) 

JO m, m„ L /=1 J 

Let us begin our discussion by defining the following summation: 
f(m x m 2 m z z)= £ exp [-Trz{am\+bm\ + cml)]- 1 , 

Z (o o o)^* s/ ^ ( fi - 3 ) 

We now make repeated use of the first identity to re-express the above function 
in the following form: 

f(m l m 2 m 3 z) = y=z £' exp [ - irazm] - irbzm\ - nm]/ cz] 

V CZ m,,m 2 ,m 3 

+ -^= J' exp [~Trbzm\--'Trmllcz] 

\l CZ m 2t m 3 

— J-_ £' exp [-nazm]- irml/cz] 

v cz m ,,m 3 

1 1 



+ 7—7= I' exp[-vm 2 x /az-<jrzbm\] 

+£' e _7rzcm 3-f-L Y' e"^ 
>A cz 1 



m 3 V m , 



+ rc-"* m Hv^v^)- 1 . (B.4) 

m 2 
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We now multiply the previous expression by z (5/2 ~° and integrate: 



=2 1 r- 5/J c- i/2 ras)f(5) 



+ 2 ,-<-i^ a -<.-iJ/a r(l(- _ _ i) 

Vc 

+- 2 =(^)-< s -^ 2 r(Ki-2))^(i-2) 

vac 

+-L f"d*2 << - W2 - n Z exp[-*rm 2 /az-irz&m 2 ] 

V AC Jo m,,m 2 =l 

+4= (""d**"-"' 2 - 0 I exp[-7razm 2 -«m|/ cz] 

V C Jo W|,m,- I 

+ -7= f°°dzz ((, - 1,/2 - ,) I exp[- ff 6zmi-7rm 2 /cz] 



+7= [°°dzz ((s - ,)/2 - 1) I exp[-7rzam 2 -7rz6mj- W m 2 /cz]. (B.5) 
All integrals can be re-expressed simply in the following form: 

j°° dz 2 a ~ x e~ z - ,/z = 2 2a+l SHt a e-**r-\a +%) J" dz (z + z 2 )- 1 ' 2 e~ w7 . (B.6) 

The right-hand side of the previous expression can be easily integrated for various 
values of a. (When a = 1, the integral yields a combination of Bessel functions 

All integrations can now be explicitly performed. When terms are collected, we 
arrive at (3.12). 

Appendix C 

In this part, we will calculate the limit of the zeta function when one of the radii 
goes to infinity. 
We start with the definition of the zeta function: 



'l?"- n iJM-< ! ?4 + w)}- , l <ci) 



508 



M. Kaku t J. Lykken / Fissioning universe 



Using the identity which flips exponentials of z into 1/z, we find: 



= 2 1 r- l/2 rG S )p^(5)+27r- ((s - ,,/2, r(Ki-l))p/3 I -^( S ) 
+2 w -« I - W2 »r(|( S -2))p/8p*- 2 f( S -2) 

+^ f'dzz^ 2 - 0 I exp[-Wf. 2 /p 2 -- I 7 2 /3 2 ik 2 /z3 




+4p f dzz ((l - 2)/2 - ,) Iexp[- } rz/j 2 /p-7rp 2 m 2 /z] 



Jo 



+MT dz /u-2)/2-i) £ exp [-,rz«7p 2 -i7p 2 ro 2 /z] 




+^ ["dzz^- 2 ^-" 1 exp[-7rzn 2 /p 2 -7rp 2 m 2 /z-ir/3 2 Jt 2 /z]. 




(C.2) 



Now we use the formula that the modified Bessei function can be written as: 
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= 2(/3pfc/ W ) ( ^ 2)/2 ^ (s „W27Ti3/cn/p) . 
Putting everything together, we find our final result (4.1). 



(C.3) 
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